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A brief introduction is given to the concept of the effective 
average action. Its dependence on the averaging or coarse 
graining scale is governed by an exact RG equation for which 
nonperturbative approximation schemes are described. This 
formalism is applied to the computation of the equation of 
state for two flavor QCD within an effective linear quark me- 
son model. Our results allow to derive the temperature and 
quark mass dependence of quantities like the chiral conden- 
sate or the pion mass. A precision estimate of the univer- 
sal critical equation of state for the three-dimensional 0(4) 
Heisenberg model is given. The exact RG formalism applied 
to the linear quark meson model is demonstrated to provide 
an explicit link between the 0(4) universal behavior near the 
critical temperature and zero quark mass on one hand and 
the physics at low temperatures and realistic current quark 
masses, i.e., the domain of validity of chiral perturbation the- 
ory on the other hand. 



I. THE CHIRAL TRANSITION AND THE 
LINEAR QUARK MESON MODEL 

The chiral phase transition is very difficult to tackle 
analytically. The main obstacles which arise are twofold. 
QCD as the microscopic theory of strong interactions is 
formulated in terms of quarks and gluons. The IR be- 
havior of strong interactions is, however, dominated by 
collective degrees of freedom like mesons. For instance, 
the relevant modes for the chirally broken phase are Gold- 
stone bosons which, in QCD, are quark-antiquark bound 
states. Furthermore, for the relevant length scales the 
running gauge coupling a s is large and a perturbative 
treatment is questionable. 

A popular way to circumvent the first difficulty is the 
use of effective field theories for the most important de- 
grees of freedom. The most prominent example is chiral 
perturbation theory based on the nonlinear sigma model 
which describes the IR behavior of QCD in terms 
of the Goldstone bosons of spontaneous chiral symmetry 
breaking. This yields a very successful effective formula- 
tion of strong interactions dynamics for momentum scales 
up to several hundred MeV or temperatures of several 
lOMeV. For somewhat higher scales additional degrees 
of freedom like the sigma meson or the light quark flavors 
will become important and should be included explic- 



itly. We will therefore rather work with an effective lin- 
ear quark meson model ||-[7|. Here the lightest mesonic 
degrees of freedom are encoded in a complex scalar field 
matrix $ which can be thought of as a (color-neutral) 

quark-antiquark composite, $ ah ~ where a, b la- 

bel the Nf light quark flavors. Spontaneous chiral sym- 
metry breaking corresponds to a scalar vacuum expecta- 
tion value ($ ab ) = a S ab with W ^ 0. For the scale- 
dependent effective action which describes the interac- 
tions of the light scalar and pseudoscalar mesons with 
quarks we make the Ansatz 

r fe [$, = J d^jz* tr [d^d*®] 

+Zy* a i$ * a + U k ($, $ f ) - \ tr ($tj + /$) (1) 

where k denotes the relevant momentum scale in a way 
defined more precisely below. The meson and quark wave 
function renormalizations, Z$ and respectively, are 
now functions of k and h is a fc-dependent Yukawa cou- 
pling. The scale-dependent effective potential is denoted 
as Uk- 

We will assume that ([!]) describes the most important 
low energy degrees of freedom of strong interactions for 
scales k below a "compositeness scale" k<$> at which the 
original gauge interaction of QCD becomes strong enough 
to trigger the formation of the light mesonic bound states. 
The effect of explicit chiral symmetry breaking due to 
non-vanishing current quark masses is represented by an 
external source j for the scalar meson field It is related 
to the average current quark mass m = (m u +TOd)/2 (we 
neglect here isospin violating effects) by 



(2) 



* Email: D.Jungnickel@thphys.uni-heidelberg.de 



where Wik& denotes the bare scalar mass parameter con- 
tained in Uk at the compositeness scale k$. The main 
approximation made in (|l|) is the neglect of higher deriva- 
tive terms as well as higher dimensional interactions 
among quarks and mesons. We note that this model is a 
generalization of the Nambu-Jona-Lasinio (NJL) model 
P,[Tof where the four-fermion interaction has been elim- 
inated in favor of an auxiliary field <E>. Beyond the NJL- 
model we allow here for a non-vanishing kinetic term for 
the scalar mesons, i.e., Z$, ^ 0, as well as an arbitrary 
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form of the effective average potential Uk for all values 
of A. 

The success of the NJL-model to describe the strong 
chiral dynamics at low energies (see, e.g., [Tot and refer- 
ences therein) serves as a further support for the use of 
(|l|). Our motivation to consider the linear quark meson 
model as an effective field theory for the IR behavior of 
QCD is, however, not only limited to the relation of (|l|) to 
NJL-models. It can be demonstrated that, similarly to 
chiral perturbation theory, also within the scalar part of 
the linear quark meson model a systematic expansion of 
observables like meson masses or decay constants in pow- 
ers of the light current quark masses is possible g. This 
allows for a successful "prediction" of the light quark 
mass ratios || and a computation of several of the Lj 
coefficients of higher dimensional operators of chiral per- 
turbation theory J?J in good agreement with the results 
of other methods. 

We will restrict ourselves here to the two lightest quark 
flavors and neglect isospin violation. Furthermore, we 
decouple the rf and the ao isotriplet. This is justified 
by their relatively large masses and can be achieved for 
Nf — 2 in a chirally invariant manner. This leaves us 
with the 0(4)-symmetric Gell-Mann-Levy model for the 
three pions and the sigma meson coupled here, however, 
to the up and the down quarks (instead of the nucle- 
ons). It should be noted, though, that our justification 
for the decoupling of the rf meson is less clear at finite 
temperature. It has been speculated Jll|,[l2j that close 
to the chiral transition temperature an effective restora- 
tion of the axial !7a(1) symmetry might take place. In 
this case the rf would become degenerate with the pions 
even in the spontaneously broken phase. In principle this 
question could be address within the linear quark meson 
model by allowing for a finite explicitly Ua(X) breaking 
scalar interaction term ~ 77(det <&+det $t) j n the effective 
average potential Uk ■ This would lead to a finite rf mass 
even in the chiral limit and the temperature dependence 
of the coupling V(T) would yield the require information 
about an approximate restoration of J7a(1) within the 
linear quark meson model. We leave this generalization 
of our model to future work. 

We furthermore note that even within this relatively 
simple effective model we have to deal with strong cou- 
plings. For instance, the renormalized Yukawa coupling 
acquires a value h ~ 6 to reproduce a realistic constituent 
quark mass of M q = hf w /2 ~ 300 MeV. Thus a nonper- 
turbative method is required which we discuss next. 



II. EFFECTIVE AVERAGE ACTION AND 
EXACT RG 

The concept of the effective average action Tk is most 
easily introduced by considering an 0(7V)-symmetric 
scalar model with real fields x a m d Euclidean dimensions 
and classical action S[x]- We define the scale dependent 
generating functional for connected Green functions as 



W k [J]=)n / V X expi-S k \x] 



d d xJ a (x) X a (x) 



(3) 

with IR cutoff scale k. Here S k [x] = S[x] + ASfcfc] with 
S[x] the classical action and 



-0^Rk(q 2 )Xa(~q)x a (q) 



(4) 



We require that the IR cutoff function Rk(q 2 ) vanishes 
rapidly for q 2 3> fc 2 whereas for q 2 <C k 2 it behaves as 
Rk(q 2 ) — k 2 . This implies that all Fourier modes x a (q) 
with momenta smaller than k acquire an effective mass 
~ k and decouple while the high momentum modes of x a 
are not affected by R k - The classical fields $ a = (x a ) — 
SWk[J]/SJ a now depend on fc, and the effective average 
action is defined as 



T fe [$] = -W k [J] + J d d xJ a {x)$ a {x) - As k m 



(•5) 



In order to define a reasonable coarse grained free en- 
ergy we have subtracted in (|^) the infrared cutoff piece. 
This guarantees that the only difference between Tk and 
r is the effective IR cutoff in the fluctuations. Further- 
more, it has the consequence that does not need to be 
convex for k > whereas a pure Legendre transform is 
always convex. (The coarse grained free energy becomes 
convex only for k — > 0.) This is important for the descrip- 
tion of phase transitions, in particular, first order ones. 
We note that the effective average action is a continuum 
implementation of the Wilson-Kadanoff block-spin ac- 
tion |E|,[l4|], however, formulated here for the generating 
functional of IP I Green functions. 

For choices of the IR cutoff function Rk like 



Rk(q 2 



Z*q 2 



e q2/k2 
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(6) 



the effective average action Tk interpolates between the 
classical action in the ultraviolet and the full effective 
action in the infrared: 



lim r fc [$] = S[$] 
iimr fc [*] =r[$] 

fc^O 



(7) 



The high momentum modes are very effectively inte- 
grated out in (0) because of the exponential decay of 
Rk for q 2 3> fc • All symmetries of the model that are 
respected by the IR cutoff AkS are automatically sym- 
metries of Tk- Hence there is no problem incorporating 
chiral fermions since a chirally invariant cutoff can be 
formulated HQ§. 

The dependence of Tk on the coarse graining scale k is 
governed by an exact renormalization group (RG) equa- 
tion fT| 
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d t T k [$] = -Tr- 



1 



-StRk 



(8) 



Here t — ln(fc/A) with some arbitrary momentum scale 
A. For the linear quark meson model we will identify A 
for convenience with the compositeness scale fc$ where 
the mesons a thought to form within QCD. The trace 
includes a momentum integration as well as a summation 
over internal indices. The second functional derivative 



(2) 



ab ( q ' q ^ 5<P a {-q)6<S> b (q') 



(9) 



denotes the exact inverse average propagator. The ex- 
act RG equation (|sj) allows to start with a microscopic 
(classical or possibly already effective) action at some UV 
scale A and successively include all quantum fluctuations 
by lowering k to zero where subsequently all physics is 
extracted. 

Eq. (|^) is a partial functional differential equation 
which, in general, will be difficult to solve. Its usefulness 
depends crucially on the existence of systematic approx- 
imation schemes. Besides conventional perturbative or 
large-A^ expansions there are two main directions which 
have been explored in the literature: an expansion of the 
effective Lagrangian in powers of derivatives 

r fe [$] - J d d x{u k { P ) + ^z k (<s>)d^ a d^ a 

+ l Yk (p)d^p + 0(d 4 )} (10) 
where p = |$ a $ a , or one in powers of the fields 

n=0 ' J \j=0 } 



X if >(*!,. 



(ii) 



(If one chooses $o as the fc-dependent VEV of the se- 
ries (O) starts effectively at n = 2.) The flow equations 

— (n) 

for the 1PI n-point functions T k are obtained by func- 
tional differentiation of cq. (§1) . The formation of mesonic 
bound states J18| , which typically appear as poles in the 
(Minkowskian) four-quark Green function, is most effi- 
ciently described by an expansion like ( pT|) . On the other 
hand, a parameterization of as in (JIG) seems partic- 
ularly suited for the study of phase transitions and will 
therefore be our choice here. 

Our Ansatz for the effective average action of the lin- 
ear quark meson model is provided by ([!]). Inserting it 
into the exact RG equation this complicated func- 
tional differential equation is converted into a set of cou- 
pled ordinary and partial differential equations for the 
fc-dependence of the two wave function renormalizations 
Z$ , , the Yukawa coupling h and the effective average 
potential U k . For the sake of brevity we will not display 
here the full set of flow equations and refer the interested 



reader to |||J. It is, however, illustrative to consider the 
flow equation for the effective average potential as an ex- 
ample. With the abbreviations u(t,p) = k~ d U k (p) where 
p = k 2 ~ d Z<f,p it is given by 
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-du + (d - 2 + 77$) pv! + 2v d ^3l^(u'; 77$) 
+l d (u' + 2pu"- m )-2i +1 Nji F)d ( 1 -ph 2 - m )} . (12) 



Here v^ 1 = 2 d+1 ir d / 2 Y{d/2) and primes denote deriva- 
tives with respect to p. We will always use in the follow- 
ing N c — 3 for the number of quark colors and d = 4. 



Eq. (12) is a partial differential equation which governs 
the flow of the effective average potential u(t, p) with 
k starting from a fixed UV scale to k = 0. Supple- 
mented by similar (but ordinary) differential equations 
for the flow of the remaining parameters 77$ = — ^ In Z$, 
rjy = —-jrlnZy and h 2 = Z^Z^ 2 h it can be solved 
numerically on a computer. 

An important nonperturbative ingredient appearing in 
all these flow equations and, in particular, also in ( |l2| ) are 
so called mass threshold functions. Full definitions of all 
threshold functions relevant for the truncation (|l|) of the 
linear quark meson model can be found in Q . A typical 
example appearing in ( [l2|) is 



M 2 



A d t (Z^R k (q 2 )) 
(2tt) 4 [P( q 2) + M 2 ] n+1 

(13) 



with P(q 2 ) = q 2 + Z^ 1 R k (q 2 ). These functions decrease 
monotonically with M/k and decay ~ (k 2 /M 2 ) n+1 for 
fc<M where M can be seen from (|1J) to be a renormal- 
ized scalar mass of the model. This implies that the main 
effect of the threshold functions is to cut off quantum 
fluctuations of particles with masses M 2 ^> k 2 . These 
functions therefore automatically and smoothly decouple 
massive modes from the evolution of the system when- 
ever k becomes smaller than the (fc-dependent) mass of 
a particle. The relevant masses typically receive contri- 
butions from arbitrarily high powers of the coupling con- 
stants of the theory and are therefore nonperturbative in 
nature. In fact, this nonperturbative dependence of the 
flow equations on the coupling constants of the model is 

already manifest in (0). The second functional deriva- 

(2) 

tive rjj. ' generically depends on all couplings present in 
the effective average action which therefore appear in 
denominators of the beta functions for these very same 
couplings. This nonperturbative dependence of the beta 
functions on the couplings is an important difference to 
Polchinski's exact RG equation p9| . It corresponds to 
a resummation of the contributions from infinitely many 
Feynman diagrams already at the level of the beta func- 
tions. It is the main reason that nonperturbative prob- 
lems like the computation of critical exponents (see be- 
low) can be carried out in a straightforward way. 
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III. INITIAL CONDITIONS FOR THE RG FLOW 

We will assume that the linear quark meson model is a 
reasonable effective description of the chiral IR dynam- 
ics of QCD for scales k below a "compositeness scale" 
/c$ where the mesons are assumed to form due to strong 
gluonic interactions. For a solution of our coupled set of 
RG equations for Z$, Zq,, h and Uk we therefore need 
initial conditions at the scale A = fc$ . We motivate these 
conditions by exploiting the relation of the linear quark 
meson model to the NJL-model. There the scalar me- 
son field $ is introduced as an auxiliary field for k — fc$ . 
This implies that its wave function renormalization Z$, 
vanishes at the compositeness scale. Furthermore, in the 
NJL-limit the effective average potential Uk is a purely 
quadratic function of $ with a positive mass squared co- 
efficient m 2 ^ ~ 1/G where G is the four-fermi coupling 
of the NJL-model. We will be somewhat more general 
here and allow for a non-vanishing by small scalar wave 
function renormalization Z$(fc$) <C 1 at the composite- 
ness scale. Moreover, we will assume that the effective 
average potential is almost quadratic for k = fc$. More 
precisely, we take a starting function Uk* which 
has its absolute minimum at the origin $ = $t = but 
is otherwise unrestricted. The two remaining initial con- 
ditions can be freely chosen as Z-q,(k$) = h(k<s>) = 1. The 
condition on the effective average potential implies that 
the evolution of the system starts at k$, in the chirally 
symmetric regime. Clearly, because of quark fluctuations 
with a large renormalized Yukawa coupling, the scalar 
mass squared term will decrease with k and finally turn 
negative, therefore spontaneously breaking chiral sym- 
metry. This typically happens at scales k x sB — 400 MeV 
well above Aqcd- This may serve as a qualitative justi- 
fication for the neglect of confinement effects in this as 
well as in the limiting NJL-model: The only degrees of 
freedom in the model which experience confinement di- 
rectly are the two light quark flavors. Yet, they acquire a 
constituent mass of the same order as k x sB before k gets 
too close to Aqcd- They will therefore decouple from 
the evolution of the mesonic system before confinement 
effects start to modify the quark propagator significantly. 

The smallness of Z$(k$) implies that all mesons have 
large renormalized masses near fc$ and therefore decou- 
ple effectively. The model is thus dominated by quark 
fluctuations for the beginning of its A; -evolution. As a 
consequence, the flow equations exhibit an approximate 
partial IR fixed point not far below fc$ More pre- 

cisely, except for the meson mass term all scalar self in- 
teractions contained in Uk divided by appropriate powers 
of the renormalized Yukawa coupling rapidly approach 
fixed points and almost all information about the initial 
conditions at fc$ is washed out. This implies an enor- 
mous predictive power of the exact RG analysis of the 
linear quark meson model which goes far beyond the chi- 
ral Ward identities and has been successfully tested |||| . 
There are only three unknown parameters left in our 



model which have to be fixed using phenomenological 
input: the (bare) scalar mass m*^ at k§ as the only 
relevant initial condition, the compositeness scale fc$ it- 
self and the average light current quark mass to. Using 
fir = 93 MeV, M w = 135 MeV and a constituent mass 
M q = 300 MeV as input we find e.g., fc$ ~ 600 MeV. 

IV. CHIRAL PHASE TRANSITION AT 
NON- VANISHING TEMPERATURE 

The introduction of finite temperature in the exact RG 
formalism is straightforward. Since we are working in 4c? 
Euclidean space time, at finite T the id momentum in- 
tegral implicitly contained in the trace in (||) is replaced 
in the standard way by a 3d integration times a Matsub- 
ara sum. Consequently, all scalar and fermion degrees of 
freedom are substituted by infinite towers of Matsubara 
modes with increasing T-dependent masses. Technically, 
one can treat T as an external parameter and evolve the 
effective average action ([j]) for each value of T from k$ to 
k = using (gj). In this context it is important that the 
initial conditions at for the RG flow discussed above 
are practically insensitive to T-effects up to temperatures 
of approximately 170 MeV S. It is crucial that the de- 
coupling of massive modes through threshold functions 
also works for Matsubara masses. For instance, eq. ( p"3| ) 
is now replaced by 

f d*q dtjZ-'Rkiq 2 )) 
J {2ttY [P(g2) + M 2 ]" +1 1 ' 

where 

q 2 = (21ttT) 2 + q 2 . (15) 

Using P(q 2 ) = q 2 + Z^Rk{q 2 ) we see that the squared 
scalar mass M 2 is effectively replaced in (|lj) by the 
Matsubara masses M 2 + {21ttT) 2 , I e ZZ. The thresh- 
old functions therefore automatically decouple all mas- 
sive Matsubara modes as k is lowered and the model 
is dimcnsionally reduced in a smooth way. This proce- 
dure yields, for instance, the full effective potential as 
a function of <f>, T and the external source j and there- 
fore the equation of state dU(($>) , T)/d$> = j which con- 
tains the T-dependence of various observables. In fig. |l| 
we show the chiral condensate ^WW) as a function of 
T/T c . Lines (a), (b), (c), (d) correspond to TO 7r (T = 0) = 
0, 45 MeV, 135 MeV, 230 MeV, respectively. For each pair 
of curves the lower one represents the full T-dependence 
of whereas the upper one shows for comparison 

the universal scaling form of the equation of state for 
the 0(4) Heisenberg model. We see a sharp second or- 
der phase transition in the chiral limit (this claim will be 
substantiated further below by demonstrating the appro- 
priate critical behavior near T c ). 
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FIG. 1. The chiral condensate as a function of T/T c . 

For non-vanishing quark masses this changes into a 
smooth crossover which becomes smoother for increasing 
Til. The values for the respective (pseudo)critical temper- 
atures are given in table Q. It is remarkable that T c for 
the largest pion mass, computed within this very simple 
effective model, is quite close to the result of current full 
two flavor QCD Monte Carlo simulations performed for 
comparable quark masses J2(J . We consider this as an in- 
dication that our model yields indeed even quantitatively 
a reasonable description of the chiral transition. This 
hold for non-universal as well as for universal properties 
of the transition as will be demonstrated below. This 
claim is further substantiated by the following considera- 
tion: Because of the relatively large sigma mass |||| , fol- 
low temperatures the linear quark meson model is well 
approximated by the nonlinear sigma model. In fact, 
our (^^E') (T) curve rather accurately agrees with that 
of chiral perturbation theory for T 60 MeV as shown 
in figure ^|. Here the solid line displays the exact RG 
result for vanishing average current quark mass rh = 
whereas the dashed line shows the corresponding result 
of three-loop chiral perturbation theory |2l],|2^] On the 
other hand, in the vicinity of T c the universal behavior of 
the model sets in which is independent of the details of 
the effective action used in our approach and also accu- 
rately described by our method as demonstrated below. 



TABLE I. The table shows the critical and "pseudocritical" 
temperatures for various values of the zero temperature pion mass. 
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FIG. 2. The chiral condensate {^'I'*! 1 ^ as a function of T//4 in 

comparison with the result of chiral perturbation theory. Here 
denotes the pion decay constant for m = 0. 

In order to demonstrate our ability to compute the 
complete temperature dependent effective meson poten- 
tial we plot in figure || dU(T)/d(f)R as a function of the 
renormalized field variable 4>r = (Z$p/2) 1//2 for different 
values of T. The first curve on the left corresponds to 
T = 175 MeV. The successive curves to the right differ 
in temperature by AT = 10 MeV down to T = 5 MeV. 
One nicely observes the convexity of the potential even 
deep in the spontaneously broken phase, i.e. for small 
enough temperatures. This property of U is rather diffi- 
cult to reproduce correctly in perturbation theory. 




<WMeV 

FIG. 3. The plot shows the derivative of the meson potential 
U(T) with respect to the renormalized field <f>n = for 
different values of T. 



The effective potential U(4>r, T) also contains the nec- 
essary information to compute the renormalized masses 



5 



(or rather inverse spatial correlation lengths) of the pi- 
ons and the sigma meson. These are obtained from U as 
second derivatives with respect to the appropriate scalar 
field components. Our results for the T-dependence of 
these two quantities are plotted in figures ^ and ^, re- 
spectively, for three different values of M^(T = 0). 
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FIG. 4. iriTT as a function of T. 
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FIG. 5. m a as a function of T. 

The solid lines in both cases correspond to the realis- 
tic value m = TOphys whereas the dotted line represents 
the situation without explicit chiral symmetry breaking, 
m = 0. The intermediate dashed lines are computed 
for m = m p h ys /10. We note that, even for the case of 
the physical current quark mass, is to a good ap- 
proximation a monotonically increasing function of T for 
all three quark masses. This implies, in particular, that 
for two light quark flavors there is no long pion corre- 



lation length in thermal equilibrium. This is interesting 
in view of recent speculations (2^j2^] that such long cor- 
relation lengths might lead to spectacular experimental 
signatures in heavy ion collision experiments due to the 
formation of large regimes of a disoriented chiral con- 
densate. A further interesting result is the dip in the 
T-dependence of the sigma mass. For vanishing tem- 
perature the sigma decays predominantly into two pions 
making it extremely wide through the large scalar self- 
interaction. From figure |s| we see that for increasing tem- 
perature the phase space for this decay decreases until 
at temperatures around ~ 100 MeV this decay becomes 
kinematically impossible. One would therefore expect to 
see a strong T-dependence of the sigma decay widths at 
least in thermal equilibrium. 

V. UNIVERSAL CRITICAL BEHAVIOR 

In addition to the non-universal results discussed 
above also universal properties of the second order phase 
transition are described accurately by the exact RG 
method. Using the equation of state allows to substan- 
tiate the claim that the chiral transition in the chiral 
limit in = is indeed of second order and to compute 
the associated critical exponents. Table [n] shows our re- 
sults which correspond to the critical exponents of the 
three-dimensional 0(4)-Heisenberg model. Our results 
are denoted by "ERG" whereas "MC" labels the expo- 
nents obtained by lattice simulations [ p5| . The agreement 
is within a few percent except for the anomalous dimen- 
sion. The latter deviation can be understood at least 
qualitatively as a consequence of the rather crude ap- 
proximation of the momentum dependence of the scalar 
propagator in (|lj). In fact, if the first non-trivial or- 
der of the derivative expansion in the scalar sector of the 
model is completed by including in the Ansatz (Q) a term 
~ y$(/o)9 M p9 M p the values of the critical exponents for 
general 0(iV)-models can be seen to improve pr| . 

Yet, the exact RG method is even capable of com- 
puting the full Widom scaling form of the equation of 
state. In fig. ^ a comparison of our results, denoted 
by "ERG", with results of other methods for the scal- 
ing function of the three-dimensional 0(4) Heisenberg 
model is shown (for notation see ||). The constants B 
and D specify non-universal amplitudes of the model. 
The curve labeled by "MC" represents a fit to lattice 
Monte Carlo data ^7j. The second order epsilon expan- 
sion [^8| and mean field results are denoted by "e" and 
"MF" , respectively. Apart from our results the curves 
are taken from p7fl . 

TABLE II. Critical exponents of the 3d 0(4)-model 
v 7 5 p T] 

ERG 0.787 1.548 4.80 0.407 0.0344 = 
MC 0.7479(90) 1.477(18) 4.851(22) 0.3836(46) 0.0254(38) 
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FIG. 6. The critical equation of state for the 3d 0(4)-Hei- 
senberg model. 

We should emphasize, perhaps, that these universal re- 
sults by no means prove that two flavor QCD is in the 
universality class of the 3d 0(4)-symmetric Heisenberg 
model. As we pointed out before there is the possibil- 
ity that the axial Ua(1) symmetry might be effectively 
restored at high temperatures in which case one would 
rather expect a (7i(2) x c7R(2)-symmetric behavior near 
the chiral transition. Our results rather confirm that the 
exact RG method is capable of computing highly nonper- 
turbative quantities like critical exponents or the scaling 
form of the equation of state with good accuracy. Assum- 
ing that the axial E/a(1) symmetry remains sufficiently 
broken even at high T we would then expect our results 
to be in good qualitative and quantitative agreement with 
the universal as well as the non-universal properties of 
the chiral phase transition of full two flavor QCD. This is 
based on the fact that our method provides a smooth link 
between two temperature regimes where the properties of 
strongly interacting matter are well understood: 

• low temperatures where chiral perturbation theory 
is expected to yield a sound description of the dy- 
namics of the Goldstone modes 

• temperature in the vicinity of T c where the univer- 
sal critical behavior of the model sets in and the 
system becomes independent of almost all the de- 
tails of the underlying microscopic description. 

The plots [I] and demonstrate that the intermediate 
range of temperature, where in principle sizeable quanti- 
tative deviations from full two flavor QCD are possible, 
is not very large. 



ACKNOWLEDGMENTS 



Applications for having provided a most stimulating con- 
ference environment. I would also like to express my 
gratitude to J. Berges and C. Wetterich for their collab- 
oration on the topics presented here. 



[1 

[2. 
[3; 
[4 
[5 

[7; 

[8 

[9 

[10 

[11 
[12 
[13 
[14 
[15 
[16 

[17 

[is; 

[19 
[20 

[21 

[22' 
[23 

[24 

[25 
[26 

[27 

[28 



S. Weinberg, Trans. New York Acad. Set. 38 (1977) 
185. 

J. Gasser and H. Leutwyler, Phys. Rep. C87 (1982) 
77. 

J. Gasser and H. Leutwyler, Nucl. Phys. B250 
(1985) 465. 

D.-U. Jungnickel and C. Wetterich, Phys. Rev. D53 
(1996) 5142. 

D.-U. Jungnickel and C. Wetterich, Eur. Phys. J. 



CI (1998) 669, flhep-ph/9606483; ). 

D.-U. Jungnickel and C. Wetterich, Phys. Lett. 

389B (1996) 600. 

D.-U. Jungnickel and C. Wetterich, Eur. Phys. J. 



C2 (1998) 557, ( |hep-ph/9704345| ). 
J. Berges, D.-U. Jungnickel, and C. Wetterich, 



HD-THEP-97-20 ( |hep-ph/9705474D . 

Y. Nambu and G. Jona-Lasinio, Phys. Rev. 122 (1961) 

345. 

J. Blinens, Phys. Rep. 265 (1996) 369. 

Pisarski and F. Wilczek, Phys. Rev. D29 (1984). 
Shuryak, Comm. Nucl. Part. Phys. 21 (1994) 235. 
G. Wilson, Phys. Rev. B4 (1971) 3174. 
Kadanoff, Physica 2 (1966) 263. 
Wetterich, Z. Phys. C48 (1990) 693. 
Bornholdt and C. Wetterich, Z. Phys. C58 

(1993) 585. 

C. Wetterich, Phys. Lett. 301B (1993) 90. 

U. Ellwanger and C. Wetterich, Nucl. Phys. B423 

(1994) 137. 

J. Polchinski, Nucl. Phys. B231 (1984) 269. 

See for instance the pro ceeding; of Lattice 97, Nucl. Phys. 

Proc. Suppl. 63 (1998) ( |hep-lat/9801024| ). 

J. Gasser and H. Leutwyler, Phys. Lett. 184B (1987) 

83. 

H. Leutwyler, Nucl. Phys. Proc. Suppl. B4 (1988) 248. 

K. Rajagopal and F. Wilczek, Nucl. Phys. B399 

(1993) 395; Nucl. Phys. B404 (1993) 577. 

K. Rajagopal, in Quark- Glu on Plasma 2, edi ted by R. 

Hwa, World Scientific (1995) ( |hep-ph/950431C| ). 

K. Kanaya and S. Kaya, Phys. Rev. D51 (1995) 2404. 

N. Tetradis and C. Wetterich, Nucl. Phys. B422 

[FS] (1994) 541. 

D. Toussa int, Phys. Rev. D55 (1997) 362, (Ihep 



lat/9607084) 



E. Brezin, D. Wallace, 
B7 (1973) 232. 



and K. Wilson, Phys. Rev. 



It is a pleasure to thank the organizers of the 5th Inter- 
national Workshop on Thermal Field Theories and Their 



7 



